Analytical solution of the Gross-Neveu model at finite density 
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Recent numerical calculations have shown that the ground state of the Gross-Neveu model at 
finite density is a crystal. Guided by these results, we can now present the analytical solution to 
this problem in terms of elliptic functions. The scalar potential is the superpotential of the non- 
relativistic Lame Hamiltonian. This model can also serve as analytically solvable toy model for a 
relativistic superconductor in the Larkin-Ovchinnikov-Fulde-Ferrell phase. 



In this paper we reconsider the simplest variant of the 
Gross-Neveu (GN) model, a 1+1 dimensional relativistic 
field theory with N species of fermions interacting via a 
quartic self-interaction 



C 



(1) 



In a previous work 0, we have found that the widely 
accepted phase diagram of this model in the large N 
limit 0] needed some revision. The dynamically gener- 
ated scalar mean field becomes inhomogeneous in a cer- 
tain region of temperature and chemical potential, a fact 
which had been overlooked so far. The four-fermion in- 
teraction then does not merely lead to mass generation 
but to formation of a kink-antikink crystal. This in turn 
reflects the presence of bound baryons in the GN model 
as can be most clearly seen in the low-density limit. The 
approach used in 0] was a numerical implementation of 
the Dirac-Hartree-Fock method (equivalent to the saddle 
point method in the functional integral approach). Al- 
though one can in principle carry out such calculations to 
any desired accuracy, it remains somewhat embarrassing 
to have to rely on numerics when dealing with a suppos- 
edly exactly solvable model. 

In the meantime we have been able to overcome this 
deficiency. In the present work, we focus on the T = 
case and explain how to construct the crystal ground 
state at any density in closed, analytical form. We expect 
that the finite temperature calculation can be done along 
similar lines but have to leave this for future work. Since 
the results of Ref. Q are fully confirmed by our new 
method we refer to this paper for more details, figures 
and a discussion of the underlying physics. For a more 
general introduction into the field of 1+1 dimensional toy 
models for hot and dense matter, see the review article 
Q . Here we concentrate on the technical details of the 
analytical solution. 



We start from the Hartree-Fock Dirac equation, 

(7 5 - i ^+7°S(x)\^(x)=^(x), (2) 



choosing the 7-matrices as follows, 



7° = -01 , 7 1 = icr 3 , 7 5 = 7°7 1 



-02 



(3) 



In terms of the upper and lower spinor components <p± 
the Dirac equation consists of two coupled equations 



dx ^ 
d 



dx 



S 



which can be decoupled by squaring, 



d 2 dS 2 
ox z ox 



(4) 



(5) 



Note that Eqs. (|4I5(I fall precisely into the pattern of su- 
persymmetric (SUSY) quantum mechanics. Let us now 
make an ansatz for S(x) based on the superpotential of 
the well-known Lame potential 0, Q 

S{x) = An = AS {Ax) . (6) 

an(Ax\K z ) 

Here, three types of Jacobi elliptic functions with modu- 
lus k appear [?J . Denoting the rescaled space coordinate 
Ax by £ , Eq. (JSJ becomes 



d 2 dS ~ 2 



The spatial period of <S(£) is 

£ = 2K 



A*' 



(7) 



(8) 
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where K is the complete elliptic integral of the first kind, 
K(k 2 ) 0. We shall choose the parameter A in such a 



2 



way that the original potential S(x) has the period a 
determined by the mean density |2j, 



7T 



1 

P Pf 



hence 



.4 



2P/K 



(9) 



(10) 



The resulting potential still has one free parameter, k, 
which determines both its shape and its size; the period 
is now fixed by the mean density. Equation JJJ) for <f) + 
can then be converted into 



d 2 

— + 2 K 2 sn 2 (£| K 2 



with 



c a 2 



(11) 



(12) 



which is recognized as the simplest case of the Lame equa- 
tion p|. The corresponding equation for differs only 
by a translation of the potential through half a period 
and thus yields an identical spectrum |6j . Our particular 
ansatz for the scalar potential was of course designed in 
such a way as to map the Dirac equation onto a soluble 
Schrodinger equation with a periodic potential. We can 
now simply use all the well-known results for the Lame 
potential. For our purpose we found Ref. .8] particularly 
useful. 

In order to determine the yet unknown parameter k we 
shall minimize the ground state energy density. In the 
Hartree-Fock approach, this is usually decomposed as 

r A/2 



(13) 

E\ is the sum over single particle energies over all filled 
negative energy states regularized by a cutoff. (For sim- 
plicity, we consider antimatter by leaving the valence 
band in the Dirac sea unoccupied Q-) E 2 is the stan- 
dard correction term for double counting of the interac- 
tion energy. Consider E\ first, transforming the Bloch 
momenta p and single particle energies u> to the corre- 
sponding quantities from the Lame equation, 



E l = -2 — 
a 1 



(14) 



where 

7r aA 
It is actually more convenient to integrate over £ , using 



dk 



e/k + k 2 -e 



d£ 2y/(l-£)(£ ~k 2 ){1 + k 2 -£) 



(16) 



Here, E is the complete elliptic integral of the second 
kind, E(k 2 ). We thus have to evaluate 



E x = -2— 





dfc 




d£ 



^£-k 2 



where now the lower limit is the band edge, 

£min 1 K 5 



(17) 



(18) 



whereas the upper limit can be inferred from Eq. i|ltj|) to 
be 



k 2 

ft. _ 



E 

1 K 



(19) 



with fc max given in Eq. I|15l) . It is necessary to keep the 
sub-leading term here, since the integral over d£ is lin- 
early divergent and the divergent part will be subtracted. 
Performing the integration in Eq. Ijl7|l yields 



E ^ 



AirKa 2 



(4E + (k 2 - 2)K) 



.^(2E + (^-2)K)l„(^ 



(20) 



The term ~ A 2 can be eliminated by subtracting the 
energy of the trivial vacuum. Now consider the double 
counting correction E2, Eq. I|13(l . in the form 



E 2 = 



I 



2Ng 2 a 2 J 



(21) 



Inserting S from Eq. JBJ and transforming to the inte- 
gration variable s — sn^, the integration can be carried 
out as follows, 



E e .s. =-2 ( % + — -L- / dxS 2 (x) ^E 1+ E 2 . j l 
Jp f 2tt 2Ng 2 aJ 



de-^^ = 2 

dn 2 ^ 



s 2 VT 



(1 - K 2 S 2 )3/2 



= -^(2E+( K 2 -2)K) . (22) 



The coupling constant is related to the cutoff via the 
(vacuum) gap equation 0, E| which reads (in units where 
the vacuum fermion mass is 1) 



Ng =InA- 
Combining Eqs. (|2T|l - lt23|) . we find 



(2E+(k 2 -2)K) In A. 



(23) 



(24) 



Upon adding E\ and E 2 and recalling Eq. JSJl, the log- 
arithmically divergent terms cancel and we obtain the 
finite, renormalized ground state energy density, 



E u _ 



4E- 



2p}K 



2E 



2)K) 

2 - 2)K) In 



(25) 



2p f nK / 
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Let us minimize this expression with respect to the mod- 
ulus k, our variational parameter. This yields the simple 
condition 



a 
1 



2p f K' 



(26) 



a transcendental equation for k. Eliminating pf from 
E lcn with the help of this relation, we finally get the 
following parametric representation of the ground state 
energy as a function of density (parameter k), 



El 

7T 



1 

47T 
1 



1 



E 
K 



2kK 



(27) 
(28) 



We also give the scalar potential S(x) corresponding to 
the optimal value of the modulus k, 



S(x) 





- 2 ) on (f 




dn h 


K 2 ) 



(29) 



S(x) interpolates smoothly between widely spaced kinks 
and antikinks (~ ±tanhx) at k — ► 1 (low-density limit) 
and the function -sin (—J for k — > (high-density 
limit). Unfortunately, it does not seem possible to ex- 
press E ren or S(x) directly in terms of pf since the rela- 
tion between pf and k, Eq. I|28|l . cannot be analytically 
inverted. 

Let us pause here and summarize what has been 
achieved so far. All we have done may be viewed as a vari- 
ational calculation of the ground state of baryonic matter 
in the GN model. Our variational ansatz amounts to gen- 
erating single particle orbits from a scalar potential S (x) 
and filling all occupied negative energy levels. Guided by 
severe restrictions from analytical solvability, we choose 
the one-parameter family of scalar potentials defined in 
Eqs. © and illUfl . The result of varying the parameter k 
is given in Eqs. I|27|) - (|29|l . This in itself does not sound 
very exciting. However, if we now compare the present 
results with those of Ref. we discover that the ana- 
lytical variational solution thus obtained agrees perfectly 
with the solution of the numerical Hartree-Fock calcula- 
tion. At all densities considered in 2] , both the values of 
the ground state energy and the shape and depth of the 
scalar potentials are indistinguishable if plotted in one 
graph. In the numerical calculation, no bias about the 
shape of the potential was put in (except for periodicity), 
since the Fourier components of S(x) were used as inde- 
pendent variational parameters. We therefore conclude 
that the true ground state happens to lie on the one- 
parameter trajectory of potentials which can be dealt 
with analytically. In order to show this without invok- 
ing any numerical results, one still has to verify that the 
ground state expectation value of -tptfj is self-consistent, 
as was done for the single baryon in [jj. This is indeed 
possible, but requires more details about the rather in- 
volved Lame wave functions Q as well as some patience 



in juggling identities for elliptic functions. In order to 
keep this paper readable, we have therefore deferred the 
full analytic proof to the appendix. 

Let us now make use of the closed formulae derived 
above to illustrate certain features of the GN crystal. If 
we go to the low- or high-density limit, it becomes possi- 
ble to systematically resolve the transcendental equation 
relating pf and k, 

k « l- 8e -^ + 32(7r + P/) e- 2 ^ 

Pf^O Pf 



1 1 



~>OC Pf Ap f Mp b r 



(30) 



The non-analytic dependence of k on p/ at low density 
is a reflection of tunnelling between the widely separated 
baryon wells. For the energy as a function of density, one 
finds 



E r . 



Pf^O 



1 

47T 



2p f 8p f 



Pf 



2 

— > oo 2ir 



2 e irp 2 f 



(31) 



In the low-density limit, the three terms correspond to 
the vacuum energy density, the contribution from the 
baryon mass (~ pMs with Mb = 2/ it) and a term de- 
scribing the repulsive baryon-baryon interaction. At high 
densities, we can identify the free massless Fermi gas 
piece, the leading perturbative correction already given 
in and the next term coming from higher order effects, 
suggesting fast convergence. In view of the comparison 
with [2| it is also instructive to determine the Fourier 
coefficients S n of S(x), 



S(x) = S n c i2 ™ x/a , 



(S- n — S n ) 



(32) 



Upon using Eq. (16) on p. 912 of Gradshteyn-Ryzhik pjj 
and correcting a misprint (ir 2 should read ir on the right 
hand side), we find the following closed expression (only 
odd n's appear), 



2 Pf 



sinh(n7rK'/K) ' 



where 



K' = K(1-k 2 ). 



(33) 



(34) 



Low- and high-density limits of our previous work f| can 
now easily be confirmed, namely 



iS„ 



iSi 



2p f 



P/ -»o sinh(n7rp/) 
1 



(35) 
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and again one finds excellent agreement with the numer- 
ical results at all densities. 

Finally, we wish to point out that the GN model at 
finite density can also serve as a solvable model for a rel- 
ativistic, inhomogeneous superconductor. Along the lines 
described in Ref. ^H, one can map the GN Lagrangian 
onto a "dual" Lagrangian which has quark-quark rather 
than quark-anti-quark pairing. All one has to do is rede- 
fine quarks into anti-quarks for left-handed quarks only. 
If one works at non-zero chemical potential, a baryonic 
chemical potential /i in the GN model corresponds to 
an "axial" chemical potential (i§ in the dual BCS-type 
model. Left-handed and right-handed fermions have op- 
posite chemical potentials, hence ^ in 1+1 dimensions 
acts like a magnetic field in 3+1 dimensions. This favors 
the appearance of the Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) phase with spatially varying Cooper pair con- 
densate 0, ^3 • K is then natural to identify the kink- 
antikink crystal of the GN model with the LOFF phase 
of the dual BCS-type model. In this sense, the present 
study may also be of some use for model studies of rela- 
tivistic superconductors. 

Appendix 

Analytical proof of self-consistency of S(x) 
We would like to show that 



S(x)=-Ng 2 J2i> a {x)^ a (o 



(36) 



The 2nd term can be simplified by partial integration and 
use of the 2nd order wave equation, Eq. {7|). Due to the 
Bloch theorem the boundary terms vanish. In this way, 
one finds 



i=2\M\ 2 - e jy^\^ + \ 



(41) 



We now insert the solution </>+ taken from the literature 
in terms of Jacobi functions 



ff (g + a) 
0(0 L 



with 



2K 



if! 



ira 
2K 



nome(n) 



(42) 



(43) 



For the upper band, a = ir/. There is a 2nd solution, <fi*j_, 
which will simply be accounted for by a factor of 2 below. 
For the definitions of the various Jacobi functions, see 0. 
Using the following addition theorem for $i 

together with standard relations between different Jacobi 
functions 0, we find 



where the sum runs over all negative energy levels corre- 
ponding to the "upper band" of the Lame spectrum (we 
are again considering antimatter). First, we have to con- 
struct normalized spinor solutions of the Dirac equation 
(0) out of the known solutions of the Lame equation (|llfl • 
We write 



i>=Af 



(37) 



and choose for cf> + a solution of the 2nd order differential 
equation l|ll|) . We recall that the relations between Dirac 
variables and Lame variables are 



£ = k 2 (oj 2 + 1) , x — k£_ , p = k/n 



(38) 



where we have used condition l|26|) . Once <p + is chosen, 
follows from Eq. 10}, 



(39) 



Let us first compute the normalization factor J\f. In 
a continuum normalization, the (spatially averaged) 
fermion density is normalized to 1 for each level, 



da; i/rip 



iaai 2 - / ddi0 + i 2 + 



(40) 



with 



\4>+\ 2 =A\ 1 



A 



cn 2 (£| K 2 
cn 2 (a|/t 2 



(45) 



(46) 



Now the d£ integration in Eq. I|41|l can be performed 
(using the same variable transformation as in Eq. (1221) as 
well as the relation I = 2K) with the result 



2\Af\ 2 A 
K 2 cn 2 (a|K 2 ) 



dn 2 (a|«; 2 ) — 



(47) 



This determines the normalization factor \N\ 2 . Let us 
now consider the scalar density. In our representation of 
the Dirac matrices, Eq. |JSJ|, 



W = -|aa| 2 (0;0- + C0 + ) 

= l -^(d s + 2SM + \ 2 . 



(48) 



In the 2nd line, we have used Eq. (|39l) . Inserting the 
expression Eq. I|45|l and performing some straightforward 
calculations yields 



2A\Af\ 2 dn 2 (a\ K 2 ) B 

W = 3 21 I 2\ S 



(49) 
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With the normalization factor \AT\ 2 determined above 
and expressing S through S(x), 



1 dn 2 (a|K 2 ) 
udn 2 {a\ K 2 ) - E/K 



S(x) 



(50) 



This shows that the ^-dependence of ipip is the same for 
each orbit, only the prefactors differ. In view of the re- 
lations [see and Eq. [fT2|l ] 



dn 2 (a|K 2 



£ -k 2 . 



we get, for negative energy states, 

Iplp = —K- 



£-n 2 -E/K 



S(x) 



(51) 



(52) 



Finally we sum over all filled states. As in the calculation 
of the ground state energy, we convert the integration 



over crystal momenta into an integration over £, include 
a factor of 2 for the twofold degeneracy of the orbits and 
employ the integration limits Eqs. H18I19|I . 



2 r £m >* d£ 


dfe 




d£ 


1 f S 


< 


— ln(A)S(a:) . 

7T 



d£- 



^{\-£){l + n 2 -£) 
(53) 



We have dropped terms of order 1/A 2 and higher, but 
of course no finite terms. Inserting the relation be- 
tween coupling constant and cutoff from the gap equa- 
tion, Eq. (|23|l . then reproduces Eq. I|36|l and proves the 
self-consistency of the scalar potential l)29[l. 
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